Neutron-star mergers in scalar-tensor theories of gravity 
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Scalar-tensor theories of gravity are among the most natural phenomenological alternatives to 
General Relativity, because the gravitational interaction is mediated by a scalar degree of freedom, 
besides the gravitons. In regions of the parameter space of these theories where constraints from 
both solar system experiments and binary-pulsar observations are satisfied, we show that binaries 
of neutron stars present marked differences from General Relativity in both the late-inspiral and 
merger phases. These strong-field effects are difficult to reproduce in General Relativity, even with 
an exotic equation of state. We comment on the relevance of our results for the upcoming Advanced 
LIGO/ Virgo detectors. 

PACS numbers: 



General Relativity (GR) has passed stringent tests in 
the solar system [ij and in binary pulsars However, 
these tests involve weak gravitational fields and/or veloc- 
ities u ^ c, so the theory remains essentially untested in 
the strong-field w ~ c regime, where high-energy correc- 
tions may appear. Because strong-field regimes arc nat- 
urally defined by systems containing black holes (BHs) 
and/or neutron stars (NS's), the final stages in the evo- 
lution of compact-object binaries provide the best oppor- 
tunities to explore gravitation at extreme conditions Q . 

Astrophysical evidence for the existence of BHs, albeit 
convincing, is still not accurate enough to discriminate 
between GR and possible alternatives. Insight in the 
nature of BHs, and on whether they are the objects pre- 
dicted by GR, may be provided by future electromagnetic 
probes (see e.g. y-lg])- Gravitational waves (GWs) pro- 
vide cleaner prospects for this task, because their genera- 
tion and propagation is largely insensitive to the compli- 
cated astrophysical environment in which BHs live. GW- 
based tests of BHs and gravity theories will be possible 
with the upcoming Advanced LIGO/Virgo detectors 0- 
12| and to an exquisite degree with future detectors such 
as LISA Ii[ill[l3|-[lll. 

Tests of gravity theories with electromagnetic and GW 
probes will also be (and to some extent, already are) 
possible with NS's. Interestingly, these objects are gen- 
erally more sensitive than BHs to the presence of extra 
deg rees of freedom in the gravity theory (see Refs. (20l - 
|2j| for some representative examples). Such behavior 
has allowed for placing constraints on gravitational theo- 
ries using observations of isolated NS's (iol - fi^ or widely 
separated binary pulsars evolving under the effect of GW 
emission (25. - ,29l | . More stringent constraints may be pos- 
sible with the expected Advanced LIGO /Virgo detection 
of binary NS's and BH-NS binaries 0, El. 

In this Letter, we consider binary NS systems and fo- 



cus on strong-field effects during the late inspiral/plunge 
until merger of the two stars (after merger effects are de- 
ferred to a future work). Although the merger is only 
marginally detectable with Advanced LIGO/VIRGO in 
GR, we will show and discuss here that if one con- 
siders modifications to the gravity theory: (i) strong 
discrepancies arise which are not captured by the sim- 
pler weak-field analyses; (ii) these effects cannot be re- 
produced within GR, even with an exotic equation of 
state; (iii) observable features will be detectable with Ad- 
vanced LIGO/VIRGO, even in the late inspiral/plunge 
and merger, unlike what happens if the gravity theory 
is unmodified GR; (iv) these features may even have as- 
trophysical implications in possible models for energetic 
electromagnetic events. 

For the gravity theory, we focus here on scalar-tensor 
(ST) theories |32H34j . in which the gravitational field is 
described by the usual tensor degrees of freedom and a 
non-minimally coupled scalar field. These viable alterna- 
tive theories of gravity can also be regarded as natural 
ones, because they include an extra scalar degree of free- 
dom in the gravitational field, as expected based on string 
theory. As well, many phenomenological gravity theories 
can be either shown to be exactly equivalent to a ST the- 
ory (c.f. for instance f{R) gravity ^,ii5j) or to contain a 
gravitational scalar besides other degrees of freedom. ST 
theories are also among the most strongly constrained 
alternatives to GR, since their history dates back to the 
50-60's with Jordan [s^, Fierz [s^j, Brans and Dicke [s^ 
(whose theory is a particular ST theory). Bounds have 
been placed on these theories with solar system experi- 
ments (3, isolated NS and binary pulsars [Ulii]. 
Stricter constraints may be obtained by detecting the 
GW signal from a gravitational collapse [13] or from vi- 
brating NS's 4l|. The remaining viable parameter space 



of ST theories is nevertheless still sizable, and these the- 
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ories are therefore a rather natural choice to investigate 
strong-field deviations from GR. In fact, we will show 
that for viable ST theories and NS binaries, strong field 
effects are possible that are qualitatively different from 
GR, and which are related to the "spontaneous scalar- 
ization" of isolated NS's in ST theories, first discovered 
in Ref. [2^. 

Methodology: Wc consider a generic ST theory 
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where k = SttG (we are setting c = 1 throughout this 
Letter), R and g are respectively the Ricci scalar and 
determinant of the metric, (j) is the extra scalar field de- 
scribing the gravitational field, 1^(0) is a generic poten- 
tial, and -0 collectively describes the matter degrees of 
freedom. Jordan-Ficrz-Brans-Dicke theory corresponds 
to the particular case uj = const, while for w = 
(w — —3/2) and a suitable potential, the theory reduces 
to metric f{R) (Palatini f{R)) gravity. Theories with 
a;((/)) — —3/2 — K/(4/31og0) are equivalent (as long as 
(j) > 1) to those studied in Ref. [20[, which give signif- 
icant deviations from GR for spherical NS's ("sponta- 
neous scalarization"). 

A somewhat simpler form for ST theories is obtained 
by re-expressing the ( "Jordan- frame" ) action ([T]) into the 
so-called "Einstein-frame" action through a conformal 
transformation gjf^ = 0(?^i/, which yields 
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where V{(f) = V {(j)) / [incfp') , and the scalar field ip is 
related to by dlog^/d^s = {2k/[3 2w(0)]}^/2^ 
posing (y9 = for = 1, this can be integrated to give 
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respectively for Jordan-Fierz-Brans-Dicke and for the 
theories considered in Ref. [2^. (Note that our is 
related to the scalar field (^sdef used by Ref. [2^ via 
V = ',5DEF/V47rG.) 



In the Einstein frame the field equations are 
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arc the scalar-field and matter stress-energy tensors in 
the Einstcn frame, and Te = T^'^gj^^. The indices of 
the matter stress-energy tensor arc raised/lowered with 
the Einstein- frame metric gE, and the relation to the 
Jordan-frame stress-energy tensor is given by T^" = 
^ j^E ^ Also, the relations between the 

matter variables in the two frames are u'^ = y/cj) (from 
the normalization condition gjf^u'^u'^ = —1); p = (fy^pE 
(from pE = u'^u'^T^j^) and p = (jy^pE (from T = ((x^Te)- 
Last, in order to have the same equation of state in both 
frames, one must have po = 'P'^Po ■ This definition, to- 
gether with the current conservation in the Jordan frame 
(Vpj'' = with j'' = pow'"), gives 



^ hJe ~ 20 
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with = Pqu'^. Therefore, solving the system given by 
Eqs. (Ill), (O, (O and ^ and transforming back to the 
original Jordan frame provides a solution to the original 
problem. This is the approach we follow in this work. 
Physical Set-up: We model the NS's with a perfect fluid 
coupled to the full field equations (|4]|6l HI) to accurately 
represent the strong gravitational effects during the evo- 
lution of a binary system. Our numerical techniques for 
solving these coupled equations have been thoroughly de- 
scribed and tested previously [42-47|. The initial data 
are evolved in a cubical computational domain defined 
by G [—350, 350] km, and we employ an adaptive 
mesh refinement that tracks the two compact objects 
with cubes slightly larger than the radius of star and 
resolution Ax = 0.5 km. 

We consider an unequal-mass binary system, initially 
on a quasi-circular orbit with separation of 60 km and 
angular velocity O = 1295 rad/s, constructed with 



LORENE [48|. The stars are described by a polytropic 
equation of state {p/c^ = Kp^) with F = 2 and K = 
123G3M|/c6. We adopt a mass ratio of = 0.937, pos- 
sible for progenitors of gamma-ray bursts [49f , and choose 
individual baryon masses of {1.78, 1.9O}M0, correspond- 
ing to gravitational masses {1.58, 1.67}M0. 

For the gravity theory, we consider y(0) = and 
w(0) = —3/2 — k/ (4/3 log 0). As mentioned, these the- 
ories are equivalent to those studied in Ref. ^ as long 
as > 1 (or equivalently, ip > 0). Because the cou- 
pling uj diverges when = 1, we initially set > 1, and 
check that this condition is satisfied throughout our evo- 
lutions. This guarantees the exact equivalence between 
our setup and that of Ref. [2^. Besides the constant /3, 
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the gravity theory is also characterized by the asymptotic 
value of the scalar field far from the source [20[. Bi- 
nary pulsars measurements require p/^AirG) ^ — 4.5 (26j . 
while the Cassini experiment constrains ipo < V'caBBini = 
2(G7r)i/V[|;3|(3 + 2wo)^/2] « 1.26 x lO-^G^/Vl/?! '('with 

cjo = 4 X 10^ [H, Moreover, from /^/(47rG) 4 

to (3/{A'kG) = —4.5, the allowed value for 959 decreases 
monotonically from 09^. to 0, again due to constraints 
from binary pulsars [26[. In our simulations, we tried dif- 
ferent allowed values of tpo , and the results do not change 
significantly when /3 is fixed. 

GW extraction and backreaction: The response of a GW 
detector far from the source is encoded in the curvature 
scalars in the physical (Jordan) frame These are 

readily obtained from the Einstein frame components as 

-04 = -Rhntrn = V'f , V'S = -Rtklrhl^ = V'f + 
i>2 = -Rikik/G = (l)i>2 + ■■■ and 022 = -Rimini = 

(^ij - Z''Z''V^V^log0/2 + ...) (with ... denoting sub- 
leading terms in the distance to the detector and l,m 
being components of a null tetrad adapted to outgoing 
wavefronts). Because far from the source one expects 
ifi = ipQ + ipi/r + 0{l/r^), with ipQ = const and ipi 
a function of x'^, and because of the peeling property 
in the Einstein frame, it is clear that -02 and '03 decay 
faster than 1/r and do not produce observable effects 
on a GW detector at infinity. However, using the fact 
that log0 = ~l3ip^ = -/3((/7§ + 2ipQipi/r) + 0{l/r^), one 
easily obtains 022 ~ P^odt^i/r. As a result, the radia- 
tive degrees of freedom (decaying as 1/r and therefore 
observable by GW detectors) are 04 (carrying tensorial 
degrees of freedom) and 022 ~ P'Podt'P'iJjl (carrying a 
purely transverse, radiative scalar mode [50|). 

It should be noted, however, that for <po — the 1/r 
radiative component of 022 vanishes. This is physically 
relevant because, as mentioned, (po is constrained to small 
values and means that for viable ST theories the purely 
transverse scalar modes couple rather weakly to GW de- 
tectors. This makes their direct detection problematic; 
which at first sight seems odd as eqs. ([6|) and ([5|) imply 
that the scalar field carries energy to infinity. (This is be- 
cause the energy flux is essentially given by the integral 
of r*"" on a 2-sphere at infinity, and T*'' ~ d^(pi/r^). It 
is easy, however, to get convinced directly at the level of 
the action ^ that these fluxes are not observable with 
GW detectors in the limit Lpo 0. The detection of 
GWs is based on free- falling test masses, so to analyze 
the detector's response one needs to look at the Jordan 
frame metric g^^/<t>{^), to which the matter fields cou- 
ple [cf. eq. ([5])]. Far from the source, in suitable co- 



ordinates one has g 
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where h^^, and S(p are small perturbations, 
we have d 



If ipo 
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= 0, 
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(fi'^) Ki 1 — pSif^ , and therefore 
at linear order. This means that 



= exp(-/ 

the motion of the detector's test masses is only sensitive 
to the tensor waves /i^^ in the limit (/Jq — ^ 0. 

Although weakly coupled to GW detectors, the scalar 
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FIG. 1: The binary's separation and the dominant mode of the ?/)4 
scalar (encoding the effect of GWs) for different values of 13. 



fluxes exert a significant backreaction on the source, be- 
cause they appear at 1.5PN order, while the quadrupo- 
lar tensor waves of GR appear at 2.5PN. More pre- 
cisely, if we consider a quasicircular binary of stars with 
masses mi and m2, and scalar charges ai and a2 [with 
a; = V AT:Gtf\/mi, where the mass scale Lpi is defined, 
as above, hY(fi=jpi^ + 'pi/r + 0{l/r'^)\, the dipole scalar 
emission is l25l 1271 l39j 



£^dipolo - 



G ( Gcsmim2 



(ai 



a2) 



(10) 



Here, Gcs ~ G{1 + aia2) is the effective gravitational 
constant appearing in the Newtonian interaction between 
the two stars, i.e. the gravitational force gets modified 
due to the exchange of scalar gravitons and becomes [sij 



F 



Gcsmim2 



The quadrupole tensor emission is instead [27|, |3S 
32G / Gcsmim2^ 



E. 



quadrupole 



5c3 



(11) 



(12) 



where v = [Gcff{mi + m2)/r]^^^ is the relative velocity 
of the two stars. Therefore, the dipole scalar fluxes, al- 
though directly undetectable, are produced abundantly 
during the binary's inspiral if the charges ai and a2 
are different, and dominate over the tensor quadrupole 
fluxes, which are suppressed by (w/c)^ relative to them. 
Numerical evolutions and comparison to GR: Our simu- 
lations confirm the qualitative features described above, 
but also highlight a more intricate phenomenology. Dur- 
ing the inspiral phase, scalar fluxes are emitted if the 
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FIG. 2; The scalar field (color code) and the NS surfaces (solid 
black line) at t = {1.8, 3.1, 4.0, 5.3}ms for I3/{AitG) = -4.5. 

scalar charges of the two stars differ. These fluxes arc 
difficult to detect directly, but clearly affect the orbital 
evolution, and their effect is thus encoded in the cur- 
vature scalar ■04, i.e. in the tensor GWs (cf. bot- 
tom panel of Fig. [IJ. Specifically, in ST-theories with 
P/{AttG) < —4.2, NS binaries merge at lower frequency 
than in GR, e.g. in Fig.[T]the plunge starts already when 
the stars' centers are ^ 52 km apart, corresponding to 
an angular velocity £7 ^ 1850 rad/s (i.e. a GW frequency 
f ^ Q,/tt ^ 586 Hz, within Advanced LIGO/ Virgo's sen- 
sitivity bands), and results in the formation of rotating 
bar (whose long-lived GW signal can be seen in the lower 
panel). Remarkably, plunges starting so early cannot be 
obtained in GR, because even with exotic equations of 
state, NS radii are constrained to i? < 14 km [5l|, so the 
interaction between the two stars cannot trigger a plunge 
until a separation ^ 2R < 28 km. 

The cause of these earlier merger is not simply the 
backreaction of the scalar fluxes (|10p (absent in GR) . In 
fact, we observe that even if we attempt to maximize 
the dipole emission ([TU]) by setting up a star with the 
maximum allowed scalar charge for a given ST theory 
(ai = Qfjnax) and one with zero scalar charge (02 ~ 0), 
the scalar field grows rapidly inside the initially non- 
scalarized star, which quickly develops a charge 012 ~ ai 
when the binary becomes sufficiently close (cf. Fig. [2]). 
This shuts off the dipole flux ([T0|) . but enhances the New- 
tonian force pulling the stars together, eq. PT|) . There- 
fore the more rapid mergers in ST theories are caused 
by the combination of dissipative [eq. ((T0| ] and conser- 
vative [eq. ([TT]) ] effects. Indeed, as a qualitative test, we 
have performed a direct integration of the PN equations 
of motion of GR with the gravitational constant G re- 
placed by Goff = G(l -I- aia2) [so as to mimic eq. (im . 
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FIG. 3: Effect of an external scalar field ipo, for I3/{AttG) = —4.5. 



with cki, Q!2 ^ 0.2 — 0.4 set to values compatible with 
our simulations], and this exercise conflrms that the en- 
hanced gravitational pull is enough to produce much ear- 
lier mergers. 

The observed growth of the scalar field and charge of 
initially non-scalarized stars that get close to scalarized 
ones can be understood in simple terms. (This phe- 
nomenon is known as "induced scalarization" (20l. |39|. 
and has also been observed for boson stars in ST the- 
ory [H^].) The scalar field of scalarized stars extends well 
beyond the radius of the baryonic matter (i^j. When 
the non-scalarized star enters this scalar-field "halo" of 
the scalarized star, it grows a significant scalar charge. 
This can be seen by studying isolated NS's in ST the- 
ory [iOl, and imposing a non-zero asymptotic value ipo 
for the scalar field, in order to mimic the effect of the 
"external" scalar field produced by the halo of the other 
(scalarized) star. The effect of (po is shown in Fig. [3l 
where we used a static, spherically symmetric code to 
calculate the scalar charge of NS's function of the 
ADM mass, for a ST theory with /^/(47r) = -4.5. As 
can be seen, even modest values of ipo induce significant 
scalar charges. This is similar, energetically, to the mag- 
netization of a ferromagnetic material immersed in a suf- 
ficiently strong magnetic field [5^. Here, the external 
scalar field makes the configuration with non-zero charge 
energetically preferred over the initial non-charged one. 

Finally, the total ADM mass of the systems we consider 
is 3.35Afo, slightly larger than the masses of observed bi- 
nary pulsars (~ 2.8— 3M0), but still plausible. This value 
was mainly chosen to maximize the deviations from the 
pure-GR case, but the effects we describe in this Letter 
are expected to arise, although in a weaker fashion, also 
for binaries with lower masses, which we plan to study in 
a follow-up work. Also, NS binaries are among the most 
likely progenitors for short gamma-ray bursts (GRBs), 
so irrespective of the properties of the currently observed 
binaries, it makes sense to investigate if modifications of 
the gravity theory can affect their physics. Our findings 
show that in ST theories the GW signal accompanying 
GRBs may be weaker (because part of the energy is car- 
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ried away in scalar waves) and different from GR, which 
would make it tricky to detect, even in the presence of an 
electromagnetic counterpart. This may have important 
implications for coincident searches of GW and electro- 
magnetic signals from GRBs. More generally, an inter- 
esting consequence of our results is that for sufficiently 
massive binaries the late-time orbiting and merger phases 
as well as the after-merger object and its dynamics can be 
significantly different from GR. This can strongly affect 
energetic events possibly associated with such mergers. 
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